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Abstract: In a large class of string and string inspired models the string excitation 
and black hole pictures invoked as an explanation of trans-GZK cosmic ray events are 
equivalent. Single particle inclusive distributions are asymptotically thermal at the 
Hagedorn temperature. The hadron multiplicities are reminiscent of multiplicities in 
heavy nucleus initiated interactions. 
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Primary cosmic rays of energies exceeding the Greisen-Zatsepin-Kuzmin (GZK) 
limit pose an exciting challenge to particle physics and/or astrophysics. In essence, 
there are no known sources within about 50 Mpc capable of producing hadrons of 
energies substantially exceeding energies of 5 x lO^^eV. Nevertheless, the trans-GZK 
cosmic rays produce air showers very similar in their structure to hadron induced 
ones. For a recent review of this puzzle, the reader may consult Sigl's article [|^]. 
In order to resolve the puzzle, we conjectured that the primaries are neutrinos, 
which are able to penetrate the cosmic microwave background (CMBR) essentially 
uninhibited. In interactions with atmospheric nuclei, however, they have CMS ener- 
gies such that ^/s exceeds the characteristic scale, M^,, of a low scale string theory p. 
Consequently, interactions are unified and string excitations produce a large cross 
section, thus generating hadron-like air showers. As an alternative, Anchordoqui et 
al, ref. proposed that ultra high energy neutrinos produce "mini black holes" as 
originally suggested by Feng and Shapere and those are responsible for the trans- 
GZK cosmic ray events. The production cross section is assumed to be proportional 
to the Schwarzschild radius squared, cf. Dimopoulos and Landsberg 

In a more general context, the connection between string theories and black holes 
has been investigated over a number of years, following the original suggestion by 
Susskind 0. A careful study of the black hole scenario was carried out by Ahn et 
al. In general, those investigations fall into two different categories. Some 

authors considered the statistical mechanics of gases made of strings, ||, |10| . Others, 



notably Amati and Russo and Damour and Veneziano |T2[ performed explicit 
calculations within the framework of some specific string models and concluded that 
highly excited string states resemble (are identical to?) black holes. 

In this article we suggest that the string excitation and mini black hole pictures 
of high energy reactions are, in essence, equivalent, provided that the black holes 
are treated quantum mechanically, as opposed to the semiclassical treatment given 
in refs. 0, |^. The conjectured equivalence holds in a broad class of theories and we 
only need rather general features of them. This is fortunate: currently there is no 
really good candidate for a phenomenologically viable string model. In fact, it has 
been repeatedly suggested that string models are perhaps not directly relevant for 
physics: one should rather think about theories of a more general category, see, for 
instance, Johnson's book for a very lucid exposition. Nevertheless, for the sake 
of definiteness, we use the term "string" in what follows. Furthermore, we exploit 
the conjectured equivalence of the two pictures in order to predict single particle 
inclusive distributions in the high energy reactions induced by trans-GZK cosmic 
rays: in principle, this prediction is testable by following the development of trans- 
GZK showers induced in the atmosphere. Here we concentrate mostly on qualitative 
aspects of the subject: a more detailed, quantitative investigation is deferred to a 
future pubhcation. 

We illustrate our point on a simple example: the excitation of a high level string 
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state. The inclusive single particle distribution considered by Amati and Russo, loc. 
cit. can be dealt with similarly. The levels of a string are labeled by a positive integer 
N such that the mass of the state is given by the formula 

~ NM^ {N > 1) (1) 

In eq. (|I]) M* stands for the characteristic string scale. In general, for ^ 1 the state 
is a a highly degenerate one. Let now O be an operator which creates a string state 
of index N from an initial state | i) . (We have in mind, for instance, string excitation 
in u-quark interactions, as conjectured in ref. 0. In that case, | i) can stand for a 
quark state and O for the fermionic current corresponding to the neutrino.) 
Omitting trivial factors, the transition probability is given by: 

Poc^(i|C"^| N,a){N,a\0\i)6{E - NM,) (2) 

a 

In eq. @ a stands for the collection of labels necessary for the full specification of 
states at level A^: one has to sum over (most of) those, since the quantum numbers 
carried by | i) and O do not fully specify the substate at level A^ to which the 
transition takes place. 

We recognize that the quantity in eq. (@), 

PM = J2 \ - NA/Q (3) 

a 

is just the (microcanonical) density matrix of the final state. (In statistical mechanics, 
a smoothing of the delta function is necessary in order to define a level density; in 
the present context, such a smoothing is automatic if the finite width of resonances 
is taken into account.) The summation over a leads to the information loss discussed 



by Amati |]T4 



It is well known that a large microcanonical ensemble can be well approximated 



by a canonical one, see, e.g. Huang's textbook |jT5l. The argument given in textbooks 
runs as follows. One establishes that the fluctuations in a canonical ensemble asymp- 
totically follow a normal distribution; the width of the distribution gets narrower as 
the energy increases. In the limit of infinite energy the normalized distribution be- 
comes a 6 function. 

A more direct approach is the following. Let E be the energy as before and H 
be the Hamiltonian of the system. One can write the microcanonical density matrix 
in operator form as follows. 

Pm = HE-H) (4) 

In the present context, the energy is the invariant center of mass energy of the 
interaction described by eq. E = ^/s and H is the mass operator of the string. 
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One then uses the following complex integral representation of the S function: 

pL+ioo 

6{E-H) = — dvex-pviE- H), (5) 

where an appropriate choice of L gives a convergence factor. Next, one takes the 
trace of both sides to obtain the partition functions. (This is sufficient: by adding 
external source terms to one obtains the generating functional of the correlation 
functions.) One has then: 

expS{E)^ — / dvexp{vE) exp{Sc{v)); (6) 

27r? JL-ioo 

We used the relation between the partition function and entropy on both sides of the 
previous equation. Obviously, exp Sc{v) is the canonical partition function at inverse 
temperature v. Next, one evaluates the integral in the previous equation by means 
of the method of steepest descents. The location of the saddle point, v = (3, is given 
by the equation: 

which is the usual thermodynamic relation between temperature and energy. Using 
the preceding equation, the leading order contribution of the integral is given by 

exp{S,{/3)-E/3)=expS{E). (8) 

The quantity S'c/3 — E(3 is recognized as the Legendre transform of the canonical 
entropy. Finally, by carrying out the Gaussian integration, one obtains a contribution 
to the entropy, which, up to an irrelevant additive constant equals to —3/2lnE/M^. 
All these formulae are valid for E/M^ 3> 1. 

An operational definition of a "large" microcanonical ensemble is that terms of 
0(1/ E) in the expression of the entropy are neghgible^. 

Here we consider level densities of the form: 

d{E) = exp S{E) ~ C {E/M,)-^-^''' exp {a{E/M,Y) ; (9) 

The quantities C, 7, a, 5 depend on the specific model considered. The contribu- 
tion of the Gaussian integral to the entropy has been taken into account. (In known 
string models, 5 — 1.) The inverse temperature is: 

^ = |^._l^ + „5f^)"M- (10) 

^ dE E \mJ * ^ ' 



^Amati and Russo loc. cit. use an adaptation of the Darwin- Fowler method and discuss some 
terms of 0(1/A^) in the framework of specific string models. For N ^ 1, the asymptotic approach 
used here and the Darwin-Fowler method yield identical results. 
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Notice that for 6 = 1, the temperature is asymptotically constant and, hence, it may 
be identified with a Hagedorn temperature, Th- Notice, however that the expression 
of the inverse temperature, eq. (|T0|) gives a temperature approaching Th from above 
From a physical point of view, this is unacceptable: the Hagedorn temperature is 
supposed to be a maximal temperature. This is due to the fact that the coefficient 
of 1/E in eq. ( ]T(]| ) is negative. The problem is corrected by taking into account 
the contribution of Kaluza-Klein excitations to the entropy. Apart from an additive 
constant, the contribution is proportional to (i?/M*)", n being the number of com- 
pact dimensions. For instance, for an open superstring in the critical dimensionality, 
n = 6 and 7 = 9/2, giving —(7 + 3/2) + n = 0. Thus, the approach to the Hage- 
dorn temperature in this case is determined by higher order terms in the asymptotic 
expansion. 

The expression of the Hagedorn temperature in terms of M^, is somewhat model 
dependent: in the calculations cited above, one has, Th = M*/(a7r), with a ranging, 
approximately, between 2 and 4 for various string models. For the purposes of 
numerical estimates in this work, we adopt a = 3: this is equivalent, of course, to 
fixing the parameter a in eq. (0). Equation ( p!0D establishes the connection between 
the string and (microscopic) black hole pictures: in fact, it means that - in a sense 
- a large and highly degenerate isolated system serves as its own thermal reservoir. 
Consequently, the terms "string" and "black hole" are used interchangeably in what 
follows. 

The asymptotic estimate given here has far reaching consequences. In particular, 
the inclusive single particle distribution can be calculated using elementary statistical 
mechanics of a canonical ensemble as it was carried out in detail for specific string 



models by Amati and Russo, ref. |Tl|. We have in the rest frame of a decaying 
resonance: 

da 1 d^fc^n 

dn = — = , , — L\ (11) 

a exp(fco/T)±l ko ^ ' 

where k^ stands for the four momentum of the particle observed in the decay and, 
asymptotically, we may put T = Th- The sign ± is valid for Fermi/Bose statistics 
of the observed particle, respectively. The quantity L is a characteristic length of 
the system. Tentatively we put it equal to a compactification radius in a symmetric 
toroidal compactification scheme. In what follows, we also put k"^ = 0. This is a 
permissible approximation except perhaps in the far infrared of the observed particle 
which, however, is not very relevant for the development of observable showers. 

We now note that the expression of dn is Lorentz invariant. Consequently, one 
can replace the rest frame expression of k^ by a Lorentz invariant quantity reducing 
to ko in that frame. This is elementary. Use the invariant: 

A-=<^ (12) 
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in place of ko in eq. ([TT|) . Here, P stands for the total four momentum of the decaying 
resonance. In a scattering process, such as 



u + q — y resonance — > (k) + X, 

one can approximate v^P^ ^ ^/s. (Here, "(A;)" stands for the observed particle of four 
momentum k.) The inclusive distribution, eq. (|11|) with the replacement k^ — > K 
shows the expected properties if evaluated in the laboratory frame. In the high 
energy limit, fermions and bosons have to be treated somewhat differently, due to 
their different statistics. 

We evaluate the invariant K in the laboratory frame, where the target quark of 
an effective mass m is approximately at rest and the incident neutrino has energy 
Ey. In the limit ^ m we have: 



Here, ko,L stands for the energy of the observed particle evaluated in the laboratory 
frame; 6 is the angle between the incident neutrino and the observed particle. 

There is a forward peak {6 = 0) both for a boson and a fermion observed. For a 
boson, we have 

1 koda 



d^k 



whereas for a fermion the term (x m/Ey can be neglected in eq. ([T3| ) and one finds 
that the forward peak remains finite: 



1 kfjda 
a d^k 



L2 
2 ' 



The different behavior of observed bosons and fermions is due, of course, to their dif- 
ferent statistics. Obviously, dn is exponentially small for a finite value of the emission 
angle 9, as expected: this is just the Boltzmann limit of the thermal distribution. 

One can evaluate the total multiplicities by means of a straight forward integra- 
tion. As stated before, neglecting the rest masses of the observed particles is a fair 
approximation. With that, one has: 

^ = ^ J ^exp(WT^)-l"I^ ^' ^ ^ 

The corresponding integral for fermions is smaller by a factor of 1/2. 

It follows that the total multiplicity of a particle of any kind is proportional to 
its statistical weight, viz. 

riB = —L^ThQe, (15) 
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and ^ 

= —L'^T^gp, (16) 

respectively. In the asymptotic limit the multiplicities are independent of the incident 
energy: an observable property of this scheme. 

As far as unconfined particles (7, W, Z, leptons) are concerned, eqs. (|l5i |l^) 
contain the final answer. However, for quarks and gluons, one has to consider their 
fragmentation into observable hadrons. Instead of carrying out a detailed analysis of 
the fragmentation, we adopt here an approximate procedure, which, however, should 
give us an idea about the order of magnitude of hadrons of any flavor produced. 

We notice that popular analytical fits to parton fragmentation functions, on 



the average, give about {nu) ^ 3 hadrons produced from each parton, see for 
a typical reference. Consequently, in order to obtain the multiplicities of hadrons 
produced, we just multiply the quark and gluon multiplicities by 3. 

Below, we list the statistical weights for each kind of particle considered, taking 
three families into account for quarks and leptons. 

9q + 9q = 72 

9g = 16 
91 + 91= 24 

9w + 9w + 9z + 9-7 = 8 (17) 

In compiling the statistical weights, we assumed that in the string regime symmetries 
are unbroken; in particular, W and Z are transverse. For the time being, we did not 
consider the various superpartners which may be produced in such reactions. 

Finally, we have to determine the characteristic length, L, occurring in eq. ([TT|). 



We use the ADD formula [|1^], assuming, for the sake of simplicity a symmetric 
toroidal compactification: 

1 /M,^'/" 



In eq. ([T8|), Mp stands for the four dimensional Planck mass and n is the number 
of compactified dimensions. Previously we presented evidence that should be 



around 80 TeV or so, [0. In what follows, we use that value of for the purposes 
of numerical estimates. However, given the uncertainties at the present stage, it is 
reasonable to conjecture 50TeV^ < < lOOTeV^. 

Using these ingredients, we estimate the total number of hadrons generated in 
the initial interaction. One has the estimate: 

Nhad ~ 3 X ((7, + 1/2 {g, + g^)) I (19) 



It is instructive to compare the multiplicity obtained from eq. ( p!9D at a reasonable 
number of extra dimensions (n = 6, 7, 8) with what an incident proton or nucleus 



- 6 - 



n 


6 


7 


8 


Nhad 


1.76 X 10^ 


4.16 X 10^ 


1.12 X 10^ 



Table 1: Hadron multiplicity for various numbers 
of extra dimensions 



would produce^. In Table 1 we summarize the hadron multiplicities estimated in the 
black hole picture, using M* = 80TeV. 

The number of prompt leptons 
can be calculated analogously. 

Let us recall that according to 
available data, the trans-GZK air 
showers are "hadron-like" : the show- 
ers develop similarly to hadron in- 
duced ones at lower energies, see Sigl, loc. cit.. Let us consider a primary energy, 
El = 3x lO^^GeV, about the highest energy observed, [Q. Assuming no new physics 
between present accelerator energies and that energy, one can extrapolate the multi- 
plicity of secondary hadrons in various reactions {pp,pp, ee). The energy dependence 
of the multiplicity is similar in all those reactions, cf. the tables in the Particle Data 
Group, |21|. Lumping all the data together, one obtains a simple fit: 

AT,;, = 2.60 - 0.18/ + 1.23f, (20) 

where / = s/lGeVl Extrapolated to El = lO^GeV in a pp interaction, one gets 
Nch ~ 170. The total multiplicity, relevant for shower development, is estimated by 
multiplying this number by 3/2, since most of the secondaries are pions. 

Clearly, the hadronic multiplicities produced by the microscopic black holes are 
larger; thus, one may conclude that the event in question was initiated by a heavy 
nucleus. Adopting a simple superposition model of nuclear interactions, see e.g. 
Gaisser ||2^ one can estimate an "effective atomic number" of the reaction by taking 
the ratio of the multiplicity produced by a black hole to that of the extrapolated 
hadronic multiplicity. All such formulae are based on a number of assumptions and 
cannot be used for obtaining quantitative estimates. 

In the following Table we summarize the effective atomic numbers, A^q, for the 
number of compactified dimensions considered here. 

Given the crudeness of the approximations used in ob- 
taining the entries of the preceding Table, this is merely a 
hint to be followed up by more careful calculations. The 
analysis of future data from detectors presently under con- 
struction will, eventually, bear out the behavior of showers 
induced by trans-GZK cosmic rays. 

We conclude with a few remarks. 



n 


6 


7 


8 


^eff 


69 


16 


4 



Table 

atomic 



2: Effective 
number 



function of n 



• We demonstrated the equivalence of a model of neutrino induced string excita- 
tions and one involving (microscopic) black holes, both purporting to explain 
trans-GZK cosmic rays. The equivalence simplifies the calculation of single 
particle inclusive distributions within such a scheme. In fact, at least asymp- 

totically. the inclusive distr ibutions can be calculated in a simple way, using 

— 6,7 are theoretically motivated: superstrings and n = 11 SUGRA, respectively. We added 

n = 8 in order to keep possible future models in mind. Low numbers of extra dimensions are 

probably excluded, see poj 
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Lorentz invariance and elementary statistical mechanics. It is to be pointed out 
that the asymptotic equivalence of the "excited string" and "mini black hole" 
pictures is independent of the tree approximation used by Amati and Russo, 
loc. cit.. It is, nevertheless, pleasing that an explicit calculation bears out the 
the black body nature of a highly excited string. 

• While we concentrated on qualitative aspects of our discussion, it is likely that 
this can be followed up by a more quantitative study, involving a MC simulation 
of the string/black hole induced trans-GZK showers. 

• There remain some questions to be clarified. In particular, the leading asymp- 
totic estimates are largely model independent, but one has little knowledge 
about model dependent corrections of 0{1/E). Very often, however, leading 
terms in an asymptotic expansion yield quantitatively correct estimates^. 

The idea of the research reported here arose in discussions with Zoltan Kunszt 
during the authors' visit at ETH in Zurich: we are grateful to Professor Kunszt for 
the hospitality extended to us and a fruitful exchange of ideas. We further thank 
Luis Anchordoqui, Will Burgett and Haim Goldberg for several enjoyable 
discussions on the puzzle of trans-GZK cosmic rays. 
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